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EXPONENTIAL ANISOTROPY OF SOLAR COSMIC RAYS 

J. w. Bieber, P. A. Evenson, and M. A. Pomerantz 
Bartol Research Foundation of The Franklin Institute 
University of Delaware, Newark, Delaware 19716 

ABSTRACT 

On 16 February 1984 a flare on the sun's invisible 
disk produced a large, highly anisotropic solar 
particle event. A novel technique, in which inter- 
planetary scattering parameters are determined purely 
from the form of the particle anisotropy, is here 
applied to energetic particle data from neutron 
monitors and the ICE spacecraft. 

1. Introduction. Recent theoretical investigations (Roelof, 1969; 
Kunstmann, 1979; Earl, 1981) indicate that the cosmic ray anisotropy 
may appropriately be expressed as an exponential function of pitch 
angle: 


f = c o + c l B exp y l M l 1 


Here, f is the phase space density, cq and ci are constants, u is 
the cosine of pitch angle, B is the magnetic field magnitude, x is the 
scattering mean free path, L is the magnetic field scale length, and q 
is a parameter which, according to quasilinear theory (Jokipn, 1971), 
is related to the slope of the power spectrum of magnetic field 
fluctuations. Provided that q and the ratio x/L do not change with 
distance, (1) is an exact solution of the steady-state Boltzmann 
equation in arbitrary guiding field configurations, including the 
Parker field. Solar particle anisotropies are, of course, not 
precisely steady-state. Nonetheless, it might reasonably be expected 
that (1) constitutes a better approximation to actual anisotropies than 
the often utilized first-order anisotropy (f = co'+ ci'u) which is 
a solution of neither the steady-state nor the time-dependent Boltzmann 
equation. 

On 16 February 1984, a flare on the invisible disk of the sun 
produced an unusual solar particle event that was recorded by neutron 
monitors on Earth (Pomerantz et al., 1984) and by the MEH instrument 
aboard ICE, which at that time was located on nearly the same nominal 
Parker field line as Earth at a distance of 0.07 AU upstream. The best 
timing of the parent flare is given by the sudden onset of radio 
emission at 0858 UT (Earth-received time) from beyond the western limb 
of the sun. This highly anisotropic event provides an ideal 
observational basis for testing which of the alternate forms of pitch 
angle distribution — first-order or exponential — best describes 
energetic solar particle anisotropies. 
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— -t N e> | ltron Monitor Observations . To determine the form of solar 
particle anisotropy in this event, the fractional increases recorded by 
8 polar neutron monitors during the interval 0900-1000 UT were 
corrected to sea level pressure. An effective value of u was 
calculated for each station. The data were then fitted to a 
first-order anisotropy and to the exponential anisotropy (1) with the 
parameter q set equal to unity. Separate determination of q from the 
available neutron monitor data was not feasible for this event. 


Results of this analysis appear 
in Figure 1. It is immediately ap- 
parent that the exponential anisot- 
ropy provides a substantially better 
description of the data than the 
first-order anisotropy. The 
goodness-of-f it parameter, x^, is ® 
nearly 10 times larger for the o 

first-order fit than for the expo- « 
nential fit. Even this poor fit for ° 
the first-order case is accomplished — 
only by resort to the unphysical -£ 
artifice of assuming negative densi- <d 
ties for u less than -0.4. Changes “ 
in p appear between the two fits be- ® 
cause the best-fit symmetry axis 
differs slightly for the two 
assumptions. 

The value of x/L corresponding 
to the exponential plotted in Figure 
1 is 2.8. However, comparable 
agreement was obtained for x/L in 
the range 2-10. The value of L at 
this time, as calculated for a 
Parker magnetic field corresponding 
to the ICE-observed solar wind speed 
of 280 km/s, was 1.4 AU. Hence, at 
neutron monitor energies (> 400 MeV) 
the scattering mean free path in 
this event was 3 AU or larger. 



- 1.0 - 0.5 0 + 0.5 + 1.0 


Fig. 1. Two alternate forms of 
pitch angle distribution are fitted 
to hourly-averaged neutron monitor 
data from Thule (o), McMurdo (•), 
Alert (o), Inuvik (■), Goose Bay (A), 
Mawson (A), Apatity (♦), Oulu (0). 


3. Spa cecraft Observations . Measurements of the flux of protons with 
energies 35-145 MeV from eight equally spaced arrival directions are 
made every 97 seconds by the MEH instrument on the ICE (formerly 
ISEE-3) spacecraft (Meyer and Evenson 1978). Since the magnetic field 
so mea sured at the spacecraft (E. Smith, private communication via 
ISEE common datapool tape), a representative pitch angle for each 
viewing direction can be computed directly by averaging the magnetic 
field observations over the period of the particle flux measurement, 
inis is in contrast to the neutron monitor analysis wherein the 

?r r iSc i 2 n »i°^J h 5 f ie ]?« wa ? a free Parameter. Figure 2 shows the 
3 o-145 MeV ICE data (10 minute average) at approximately 0955 UT on 

16 February 1984, near the maximum of the event at these energies. The 
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sectored data obviously do not follow a straight line (first-order 
anisotropy) but the exponential (1), with a choice of q=1.2 and 
x/L=2.8, provides an excellent fit, as shown by the plotted curve. 



Fig. 2. Sectored ICE particle data (circles) are compared to an 
exponential anisotropy (curve) with x/L = 2.8 and q = 1.2. 


The persistent high anisotropy of the event and the long mean free 
path determined from the shape of the pitch angle distribution suggest 
that the entire time profile is a reflection of continued injection at 
or near the sun. We have modeled the time profile using coronal 
diffusion (Reid, 1964) and interplanetary focused diffusion (Bieber et 
al., 1980) and find excellent agreement. A key result of this analysis 
is that the flow at the position of the spacecraft is at all times 
approximately steady state. Changes with time are nearly all due to 
changes in the source strength at the sun. 

Such a steady state flow may be used to look for changes in the 
scattering parameters which determine the shape of the pitch angle 
distribution. A related analysis conducted by Bieber and Pomerantz 
(1985) uses long term averages of neutron monitor data. The shape of 
the pitch angle distribution may be characterized many ways, one of 
which is to consider ratios of the amplitudes of harmonic components 
determined by Fourier analysis of the counting rate as a function of 
arrival direction. Figure 3 shows the average flux, ratio of second to 
first harmonic r 2 /ri, and ratio of third to first harmonic r^ln 
as a function of time. Note the distinctly different values of the 
harmonic ratios before and after the data gap. This jump implies a 
reduction in the mean free path for scattering by almost a factor of 
two, while q remains nearly unchanged. It is likely that this change 
is related to the change in the properties of the solar wind which is 
indicated by the increase in the magnetic field magnitude also shown in 

Figure 3. 
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Fig. 3. Time profiles of average 
particle flux, harmonic ratios 
r 2 /rj and r (see text), and 
magnetic field B observed at ICE. 
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4. Conclusions . Local 
measurements of important 
parameters of interplanetary 
particle propagation can be made 
with high time resolution if 
particle anisotropies can be 
measured with sufficient 
statistical accuracy. Analysis 
of energetic particle data for 
the event of 16 February 1984 
from neutron monitors and from 
ICE leads to the following 
conclusions: (1) Particle 

anisotropies are better 
described by an exponential 
function of pitch angle than by 
the commonly assumed first-order 
anisotropy. (2) Fitting to the 
exponential anisotropy allows 
both the scattering mean free 
path x and the parameter q, 
which characterizes the 
dependence of the scattering 
rate upon pitch angle, to be 
determined. For this event, it 
is found that x is greater than 
3 AU and q is approximately 1.2. 
(3) An apparent change in the 
scattering mean free path 
occurred in connection with a 
small change in the magnitude of 
the interplanetary magnetic 
field. 
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